On the quark mass dependence of 
nucleon— nucleon S— wave scattering lengths 

Joan Soto and Jaume Tarrus 

Departament d'Estructura i Constituents de la Materia and Institut de Ciencies del Cosmos 

Universitat de Barcelona 
Diagonal, 647, E-08028 Barcelona, Catalonia, Spain. 

e-mails: joan.soto@ub.edu, tarrus@ecm.ub.es 

May 9, 2012 

Abstract 

In the framework of a Chiral effective theory with dibaryon fields, we calculate the 
pion mass dependence of the inverse scattering length of the nucleon-nucleon system 
in the 3 Si channel at order (m^/A^) x (m l /2 m% 2 /^fl) n for all n > 0. We show 
that certain sets of potentially large higher order contributions vanish. We discuss the 
difficulties of extending the proof to the 1 S'o channel. We apply our results to chiral 
extrapolations of current lattice data. 

1. Introduction 

The size of the nucleon-nucleon S-wave scattering lengths is larger than expected from 
standard arguments of chiral counting [1, 2], and understanding their values from QCD is 
still a major challenge [3]. Lattice calculations at physical light quark masses are very costly 
(although they have recently been carried out for some observables [4, 5]) and, hence, the 
use of chiral extrapolations will be needed for some time in order to obtain reliable estimates 
[6, 7). 

The quark mass dependence of low energy observables can be obtained from suitable chiral 
effective theories. The chiral effective theory for the nucleon-nucleon system was proposed 
by Weinberg in Refs. [1, 2]. The fact that the S-wave scattering lengths are unnaturally 



large, together with other problems of the original proposal related to renormalization and 
consistency with the chiral counting [8], led to the so called KSW (Kaplan, Savage, and Wise) 
approach [9, 10] (see [11, 12, 13, 14, 15, 16, 17, 18] for alternative approaches, and [19] for a 
review on the original approach). It was soon realized that the introduction of dibaryon fields 
in the effective theory was a very convenient way of implementing large scattering lengths 
[20, 21]. Nowadays, a nucleon-nucleon effective field theory (NNEFT) with dibaryon fields 
has been used to calculate the phase shifts in the 1 S and 3 Si- 3 Di channels up to next-to- 
next-to-leading order [22, 23], providing results similar to the KSW approach [24, 25], with 
more economical expressions. It was already noticed in [24] that beyond next-to-leading 
order (NLO), part of the calculation must be organized in powers of ^Jm n / A X) m w being 
the pion mass and A x a typical hadronic scale (say A x ~ m p ~ 770MeV), rather than in 
powers of m 7T /A x [23]. It is in fact an accident due to Wigner symmetry that the would-be 
O (jrh/ 2 / A x 2 ^j correction vanishes [26, 23]. In addition, it was pointed out in Ref. [27] 
that the terms giving corrections y / m vr /A x were generically large. In this paper we show 
that these terms can be summed up in the 3 S\ channel and, furthermore, that they give a 
vanishing contribution to the scattering length. This allows us to provide a reliable chiral 
extrapolation formula for the inverse scattering length including terms up to order m^ 2 /A x 2 , 
m q being the average light quark masses. Let us recall that the quark mass and the pion 
mass are related by = 2B m q , where B is a low energy constant related to the quark 
condensate. Unfortunately neither the arguments that allow the resummation nor the proof 
that the effect vanishes apply to the 1 S channel. 

The paper is organized as follows. In section 2 the NNEFT with dibaryon fields is briefly 
reviewed. In section 3 we argue that exchanges of potential pions in loops with a radiation 
pion have to be resummed. We show that they give a vanishing contribution in the 3 Si 
channel. Section 4 is devoted to obtaining expressions for the inverse scattering lengths up 
to m 3 / ' A x terms. In section 5 we compare our results with the available lattice data. We 
close with a discussion and conclusions in section 6. 

2. NNEFT with dibaryon fields 

Our starting point is the effective field theory (EFT) for the Nb = 2 (Nb being the baryon 
number) sector of QCD for energies much smaller than A x , proposed in Ref. [22] . The distinct 
feature of this EFT is that, in addition to the usual degrees of freedom for a NNEFT theory, 
namely nucleons and pions, two dibaryon fields, an isovector (Df) with quantum numbers 
1 S and an isoscalar (D v ) with quantum numbers 3 Si, are also included. Since mjv ~ A x , 
Tfi jv being the nucleon mass, a non-relativistic formulation of the nucleon fields is convenient 
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[28]. Chiral symmetry, and its breaking due to the quark masses in QCD, constrain the 
possible interactions of the nucleons and dibaryon fields with the pions. The Nb = and 
Nb = 1 sectors are the usual ones and will be needed only at leading order (LO). 

The Nb = 2 sector consists of terms with (local) two nucleon interactions, dibaryons, and 
dibaryon-nucleon interactions. The terms with two nucleon interactions can be removed by 
local field redefinitions [21] and will not be further considered. The LO terms with dibaryon 
fields and no nucleons in the rest frame of the dibaryons read 

1. 



£o( P ) = 2 Tr 



^ (-ido + ^) + .Dj (-i^, + 5^,) ^ + ic w (^Tr - ^.c.) , (2.1) 



where D s = D a s r a and S r m , i — s,v are the dibaryon residual masses, which must be much 
smaller than A x , otherwise the dibaryon should have been integrated out as the remaining 
resonances were. The covariant derivative for the scalar (iso vector) dibaryon field is defined 
as d D s = d D s + \[[u, d u], D s \, 

4^ =%iWWiV*)/\ a + %ivVVViV)Dt a+ 

A A ^ 

+ f±{N^T 2 &<j 2 N*) ■ D v + ^=(N T r 2 a 2 aN) ■ D\ , 
v 2 v 2 

1/2 

with A S ,A V ~ A~ ' . The NLO pion-dibaryon 

£ 0{p2) = Sl Tr[D s (uM ] u + v)Mu ] )D\] + s 2 Tr[Dl(uM ] u + u ] Mu r )D s } + 
+ v x 3l ■ D v Tr[u ] Mu ] + uM ] u] + ■■■ , 

where M. = m q I. The s,, % — 1,2, and V\ are low energy constants (LEC). We have only 
displayed here the terms which will eventually contribute to our calculations. The complete 
list of operators is given in Appendix B of Ref. [23]. The tree level dibaryon propagator 
expression i/(—E + 8' m . — irj) gets an important contribution to the self-energy due to the 
interaction with the nucleons as discussed in Ref. [22]: 

i = s,v, (2.4) 



-E + 6L 



p = ^ErriN, which is always parametrically larger than the energy E. The size of the 
residual mass can be extracted computing the LO amplitude using the propagator (2.4) and 
matching the result to the effective range expansion, 

C-z^-^r. * = s > v > ( 2 - 5 ) 

ira 1 A x 



where a 1 , % — s ,v, are the scattering lengths of the 1 S'o and 3 Si channels respectively. As a 
consequence — E + S' m . in the full propagator can be expanded for p ~ m n , and hence the 
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Figure 1: Example diagram of enhancement of a potential pion inside a radiation pion loop. 

X'X-JL-IL- X- XX -xxx- 

Figure 2: Potential pion exchanges in the 1 5o channel can be approximated by contact interactions and 
resummed into an effective vertex when the external momentum is bigger than the pion mass. 



LO expression for the dibaryon field propagator becomes , 
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. i = s,v. (2.6) 
Afm N p 

The expanded terms can be taken into account through an effective vertex. Moreover, Eq. 
(2.4) implies that the dibaryon field should not be integrated out unless p <C S' m ., instead of 
E <C 5' m . as the tree level expression suggests. 

In order to calculate the scattering lengths, we need the nucleon-nucleon amplitudes at 
zero energy. Following Ref. [23], we will first match NNEFT to pNNEFT, an effective theory 
for E <C m n and p < m^, and then match pNNEFT to pNNEFT, the pionless EFT with 
dibaryon fields for p <C [21], from which we can easily identify the scattering lengths. 



3. Potential pions in loops with radiation pions 

pNNEFT is obtained from NNEFT by integrating out nucleons of energy E > m v and pions. 
Among the latter there are the so called radiation pions, namely pions with q° ~ q ~ m v 
that interact with nucleons of E ~ m n and p ~ ^Jm^ra^. We discuss in this section this 
particular class of contributions to the matching calculation. 

The lowest order diagrams involving radiation pions are depicted in Fig. 6. When a so 
called potential pion, namely a pion with q° ~ m n and q ~ y/m^m^ in this case, is added to 
one of those diagrams, for instance as in Fig. 1, a parametric suppression of only ^/m^/ 'A x 
occurs [26, 23], which numerically turns out to be 0(1) [27]. It is then necessary to sum up 
these kinds of contributions. 



3.1 Loop resummation 

Let us consider the exchange of n potential pions between two nucleon lines. If we project 
it to the 1 S'o channel, the three-momenta coming from the vertices of each potential pion 
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Figure 3: Resummation of potential pions in the dibaryon-nucleon vertex. 

exchange contract between themselves. Note that this is not the case if we project to the 3 Si 
channel, where a three-momentum from one of the vertices of a given potential pion exchange 
may get contracted with a three-momentum of a neighboring potential pion exchange vertex. 
If these n-pion exchanges are in a loop with a radiation pion, then the three-momentum 
in the denominator of the potential pion propagators dominates over the pion mass and the 
pion energy. As a consequence, the potential pion exchanges collapse into a local vertices 
(contact interactions) with a coupling constant g\/(2f%), where g& is the axial pion-nucleon 
coupling constant and f n is the pion decay constant as defined in Ref. [23]. Again, this 
is not so in the 3 Si channel, where even at very large momentum transfer the potential 
remains non-local (i.e., it does not reduce to a contact interaction). In the left hand side of 
Fig. 2 we depicted the first terms in a series of diagrams with an arbitrary large number of 
potential pion exchanges. Using the previous reasoning we can collapse the potential pion 
exchanges into local vertices obtaining the diagrams on the right hand side. In dimensional 
regularization the result for the first few terms is 

where we have taken the external energy to be — q°, and a is defined as 

9A m N 

Naively we would expect each bubble to suppress the diagram by a factor of sjm^j K x . 
However a more careful analysis shows that the actual size of each bubble is in fact ^frn^j a ~ 
1.19, which is of order O (1), and hence the series should be resummed. The result of the 
resummation can be cast as an effective energy-dependent four-nucleon vertex with coupling 
constant 

^/ = «4— ?r^=- (3 - 3) 



VI a + Jq 



An analogous resummation has to be done for potential pion exchanges in the nucleon- 
dibaryon vertex of Fig. 3. Following the same procedure as before, we obtain an energy- 
dependent effective nucleon-dibaryon vertex, 

ex 

•U./f ,1, ==. (3.4) 
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Figure 4: Inside radiation pion loops the 1 5'o receives an additional self-energy contribution. 





(a) (b) 
Figure 5: Order 0(m^/A x ) contributions to the dibaryon residual mass. 

Furthermore, using the effective vertex of Eq. (3.3) we can construct the self-energy 
depicted in Fig. 4, which inside radiation pion loops turns out to be of order O (1) and 
thus has to be included in the LO propagator (2.6). The following expression for the 1 S'o 
propagator inside radiation pion loops is obtained 

&ll + -^=\. (3.5) 



Note that in order to have a 1 5'o nucleon-nucleon state in a loop with a single radiation 
pion, the initial nucleon-nucleon state must be in the Z S\ channel. This procedure can then 
be applied to the calculation of a Sl , but not to the calculation of a So . This is due to the 
fact that in the last channel the contact interaction is replaced by a non-local potential that 
turns out to be singular, and therefore cannot be straightforwardly used in a Lippmann- 
Schwinger equation; see Refs. [11, 12, 13, 14, 15, 16, 17, 18] for discussions and possible 
solutions. 




Figure 6: Order 0(m^/ A x ) contributions to the dibaryon residual mass with one radiation pion, that 
cancel due to Wigner symmetry. 
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3.2 Cancellation of the contributions to a Sl 

Making use of the new effective vertices obtained by resumming potential pion exchanges, 
two new diagrams contributing at LO to a iSl , shown in Fig. 5, are found: 

Aa = 8A 2 V (^5(1/4, 1) - a 2 B(l, 1) + aB(3/4, 1) - 5(1/2, 1)) , (3.6) 



A, 8.4;. ( £(1/2. 1) - nL>(3/ I. 1) +<\ 2 B(Q.2) - <\ :i B( 1/ 1.2) n ( 'L'(1.2) fa 7 5(5/4,2) 

(3.7) 

The definition of B(fa,fa) can be found in the Appendix. These contributions are of the 
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same order as the diagrams in Fig. 6. Those diagrams were counted as O (ml/ 2 / A 
Ref. [23]; however they are proportional to y/m^/a ~ 1. Then analogous to what we did 
previously, we should count diagrams in Fig. 6 as O (m^/A x ). The sum of these diagrams is 
known to cancel due to Wigner symmetry; however since the third one is already included 
in Fig. 5b we should add the first two to Eq. (3.7) in order to get the complete result at 
0(ml/A x ), 

As = -8^-5(5/4,0) . (3.8) 

a 

The sum of these three contributions (A a , Ab, A s ) adds up to zero, which can be checked 
by making use of the relation 

B (fa - 1, fa) = B {fa, fa-l) + a A B {fa, fa) . (3.9) 

This is at first sight a surprising result. The interaction of nucleons with potential pions 
spoils the arguments that led to the proof that the sum of the diagrams in Fig. 6 vanishes 
as a consequence of Wigner symmetry [26]. Yet, since the contact four-nucleon interaction 
we obtain is only used in the 1 S channel, it could well be replaced by a Wigner symmetric 
one with no effect in our calculation, and hence the arguments of Ref. [26] would still apply. 
Nevertheless, as it will become clear soon, the actual reason for the cancellation is that the 
contact four-nucleon interaction can be removed by the following local field redefinition of 
the dibaryon field: 

D s -+ D s - tSt^P^N, (3.10) 

where P^ s " = ( " T2 ^^?' rc ^ , is the projector to the 1 S'o partial wave. Indeed, we have checked 
that the resummation of potential pion exchanges in the diagrams of Fig. 7, in which Wigner 
symmetry is violated by the cross and bullet vertices, also vanishes. 

As we have mentioned in the previous section, the resummation cannot be carried out for 
the analogous diagrams for a ls °. However, it is likely that the perturbative expansion also 
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breaks down in this channel due to numerical factors coming from loop integrals. Hence, 
any prediction for the quark mass dependence of a~ s ° in terms of a perturbative expansion 
has to be taken with caution, because it could be missing large corrections. 

Part of the reasoning we have used in the 3 S\ channel can be adapted to discuss the 
result for the diagrams with a single potential pion exchange in a loop with a radiation 
pion in the 1 Sq channel. In this set of diagrams, the radiation pion three-momentum in the 
denominators of the loop integral can be neglected in front of any of the nucleons or potential 
pion three-momenta, so the potential pion three momenta in the pion-nucleon vertices must 
end up contracted between themselves, and hence we are left with a situation analogous to 
the one in the 3 Si channel. At this point we can approximate the potential exchange by 
a four-nucleon contact term, following the same reasoning as for the contributions to a 3Sl . 
The contact term can then be eliminated by a field redefinition analogous to Eq. (3.10) for 

— * 

the D v dibaryon field. We then conclude that the sum of the diagrams in Fig. 6 with a single 
potential pion insertion must also vanish in the 1 S channel. This result is in contradiction 
with those of Refs. [25, 23], where this class of diagrams with one potential pion inside 
a radiation pion loop was found to be non-zero. We believe that this is a consequence of 
double counting certain diagrams. In particular, the last diagram in Fig. 17 of Ref. [25] is 
already included in the first one. According to our calculations this error would lead to the 
result presented in Ref. [25, 23]. 

4. Scattering Lengths 

In this section we sketch the matching between pNNEFT and NNEFT, and between pN- 
NEFT and pNNEFT in the light of results of the previous section. In particular we focus 
our efforts on obtaining NLO expressions for the residual mass and the dibaryon-nucleon 
vertices' low energy constants. With these expressions we write the scattering lengths up 
to m 3 T /A^ c terms. Results for the 1 S channel have to be taken with caution due to possible 
large corrections from multiple potential pion exchanges in the loops with a radiation pion, 
as explained in the previous section. The following subsections are rather sketchy. We refer 
the reader to section 4 of Ref. [23] for details on the matching procedure beween NNEFT 
and pNNEFT, and to section 6 of the same reference for details on the one between pNNEFT 
and /NNEFT. 

4.1 Matching pNNEFT with NNEFT 

In the one-nucleon sector, pion loops produce a shift in the nucleon mass, Sm^, that intro- 
duces a quark mass dependence. We can reshuffle 5m^ into the dibaryon residual mass by 
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a) 



b) 



Figure 7: Order 0(m^/A^) contributions to the dibaryon residual 



mass. 



local field redefinitions. The expression for 5m N can be found [29], and up to O (m^/A^J 
contributions it reads 



5m 



N 



-4cim n 



2 3g A 



327T/2 * 



(4.1) 



In the NLO pion-dibaryon Lagrangian (2.3) the residual mass gets O (m£/A x ) contributions 
proportional to the quark mass. Additional O (m^/A^j contributions come from the dia- 
grams in Fig. 7. Adding up all the contributions we obtain the formula for the residual 

mass, 



o v i 2 . , / 9a \ m l A l . f 9a\ m\A v 
" ~B~ mn ~ ~ ~ '' ~ ~ ~ 



g A \ ml A s 



(4.2) 



PJ 8tt A v 



In the two-nucleon sector only the one-pion exchange is relevant at this order, which pro- 
duces the well known one-pion exchange potential. 



4.2 Matching pNNEFT with /NNEFT 

The next step in order to evaluate the scattering length is to build a theory valid for p < S m . 
This is achieved by integrating out the nucleon three-momenta of order m^, which leads to 
the so-called pionless nucleon-nucleon EFT. Non-local potentials can be expanded in powers 

2 

of and become local. Self energies in Fig. 8(a) can be expanded, giving contributions to 
the dibaryon residual mass as well as time-derivative terms. The latter can be reabsorbed 
by field redefinitions of the dibaryon fields. The dibaryon-nucleon vertex gets contributions 
from the diagrams in Fig. 8(b). Recall that the one pion exchange potentials in Fig. 8 
correspond to potential pions with q ~ m^ and not to potential pions with q ~ y/m n m^, 
such as the ones considered in section 3. 

The contribution to the residual mass from the first diagram in Fig. 8(a) is of order 
0(m^/A x ). This diagram contains a divergence proportional to the quark mass which is 
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+ ... 

a) b) 

Figure 8: Contributions to the matching between pNNEFT and /NNEFT. (a) Leading and NLO contribu- 
tions to the residual mass, (b) NLO correction to the dibaryon-nucleon vertex low-energy constant. Dashed 
lines represent the one-pion-exchange potential. 

renormalized by the counterterm of the same order proportional to the quark mass in Eq. 
(4.2). The second diagram in Fig. 8(a) is NLO, O (m^/A^): 



6L° = 6L +2^ml-8 Cl ml 

-DO 



2 g\ f m n m N \2 f m\ 



5 LO = 5 , + 2 Sl + 8 2 2 _ 2 _ A 29\ ( ^7U N y f m? 



(4.3) 



rNLO 



fl V 4vT ) ~ V /i 2 

c- - + (§) iff + - (i) it + * (I) 2 ra 3 *- < 2 > ■ 

(4.4) 

The dibaryon-nucleon vertex up to NLO gets only one new contribution from the first 
diagram in Fig. 8(b). Defining A 2 iNLO = A 2 A NLO ,i = s,v, 

_ g\ m n m N 

A*zo-2^-^. (4-5) 

Note that the parametric suppression of an extra pion exchange in the diagrams of Fig. 8 
is m n mN I 'A 2 X whereas the one in the diagrams in Figs. 2 and 3 is mj 2 m^ 2, j A 2 X . Hence the 
resummation of diagrams in Fig. 8 is less important than the ones in Figs. 2 and 3, a feature 
that justifies why in our power counting scheme they need not be resummed. 

Now we have all the ingredients to write the expression for the scattering lengths: 



a-. 



NLO 



A nlo ) + ^jj ,i = s( 1 S ),v( 3 S 1 ). (4.6) 



m N Af m N A 2 
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,1 



g A m N \ 2 log(2) 
/I 



3 S 1 



27r((si+s 2 )/.B -8ci) 
m N A1 

2tt(i;i/Bo-8ci) 



j «„, A /-> I a\ 



+ 



2c s 



i 



16 mAr /2 A2 T g A A s A v 



+ 



g A A s A l 



i) 



9a (si+s 2 )/B -8c 1 
4/1 Al 



1287T 2 



' \ ffj (vi/Bq-8ci) i ^^ gjmjy ^ 6+131og(2) 



A?, 



4/1 



256tt 2 



Table 1: Independent free parameters in terms of the effective theory low energy constants. 



m n (MeV) 


a So (fm) 


a Sl (fm) 


353.7 


0.63 ±0.50 


0.63 ± 0.74 


492.5 


0.65 ±0.18 


0.41 ±0.28 


593 


0.0 ±0.5 


-0.2 ± 1.3 


390 


0.11818$ 


0.052±^1 



Table 2: Lattice data point used to fit the scattering lengths. The first three data points are from Ref. [33] 
and the fourth one is from Ref. [34]. 



5. Comparison with lattice data 

The expressions for the scattering lengths can be rewritten to collect all the parameters into 
three independent ones, 

The expression obtained is quite simple and emphasizes the m n dependence. The relation 
of the ( parameters to the low energy constants of the EFT can be found in Table 1. The 
expected sizes of these parameter are, (n ~ O (m^/A x ), Q2 ~ O (1/A X ) and C«3 ~ O (l/A^,). 

The first lattice QCD calculation of the nucleon-nucleon scattering lengths was per- 
formed by Fukugita et al [30, 31] in the quenched approximation with Wilson quark action. 
More recent studies using the quenched approximation have been carried out by Aoki et al 
[32]. The NPLQCD Collaboration has performed unquenched calculations in mixed-action 
(domain wall-staggered) [33] and anisotropic clover-quark action [34]. We fitted the lattice 
data of the NPLQCD Collaboration (see Table 2). Unfortunately all data points are above 
or close to 350 MeV, a scale beyond which it is not clear that chiral extrapolations for the 
nucleon-nucleon system are still valid. Thus the obtained results have to be taken with 
caution. We forced the expressions for the scattering lengths to reproduce the experimental 
values at the physical pion mass, a So = —23.7 fm and a Sl = 5.38 fm. This allowed us to 
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LO 


Xd.o.f 


Ci(MeV) 


( 2 (MeV- 1 ) 


1 S'o 

3 s, 


3.74 
0.91 


-126 
-98 


0.67- 10- 3 
1.59 • 10-3 



Table 3: LO fit parameters results. 



solve one parameter as a function the remaining ones, we chose to solve Q\. The remaining 
parameters have been obtained by minimizing an augmented chi-square distribution [35] for 
each scattering length. The augmented chi-square distribution is defined as the sum of the 
chi-square function with a set of priors for every one of the free parameters to be fitted, 



Xaug 



2,2 

Xtf ' X 



1 



n 
3=1 



prior i 
N 



Xprior 



1 ^ (In ||zfc|| - In \\x k \ 



= -Y 



k=i 



In 2 Ri 



- l So, 3 Si , 

(5.2) 



where a) and 5 a i stand for the value of scattering length and its uncertainity at the pion 
mass m n j respectively, n is the total number of lattice data points. Furthermore, x k refers 
to the free parameters, N being their total number. The free parameters are Q 2 at LO, 
and Q 2 and at NLO. The prior information is obtained from naive dimensional analysis. 
For instance, if the parameter x k is of order 0(1), we would expect it to be in the range 
0.1 < \\xk\\ < 10, which translates to setting ln(||xfe||) = and \n(R k ) = 1 for the kth 
parameter. We have taken logarithms in the prior functions to achieve equal weights for 
the subranges 0.1 < \\x k \\ < 1 and 1 < \\x k \\ < 10. For ( i2 , priors are set to ( i2 = j- and 
In (Rq 2 ) = 1, and for ( i3 to ( i3 = and In (Rq 3 ) = 1. The plots corresponding to the fits of 
the leading and next-to-leading order expressions of the scattering lengths as a function of 
the quark mass are displayed in Fig. 9. The chi-squared distribution per degree of freedom 
is defined as 

, 2 

, t= 1 S , 3 S 1 . (5.3) 



2 I {a\ 

Xa\d.o.j n _l_ 

i=i 



/// 



flj)' 



The values obtained for the parameters and the chi-squared per degree of freedom are 
collected in Tables 3 and 4. The values obtained for ( s i and ^3, % = s,v, at NLO are 
on the limit of what we would consider natural size. This could indicate that significant 
cancellations occur at the physical pion mass in order to produce the observed values of the 
scattering lengths. Note that the fine tuning increases with the precision of the expression 
used. 

The m g -dependence of the scattering lengths has been studied previously in Refs. [6, 36] 
using numerical solutions to the Lippmann-Schwinger equation with potentials obtained 
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Figure 9: Plots of a So (left) and a Sl (right). The solid and dashed lines correspond to the LO and NLO 
respectively. The triangular dot in the a Sl figure corresponds the physical value of the scattering length. 
In the a s ° figure the physical point is out of scale. 



NLO 


~X.d.o.f 


Ci(MeV) 


UMeV- 1 ) 


UMeV- 2 ) 


'So 


2.4 


-246 


4.56 • 10- 3 


9.21 • 10- 6 


3 s 1 


0.4 


-155 


3.83 • 10-3 


10.1 • 10~ 6 



Table 4: NLO fit parameters results. 



from Weinbergs's power counting, and in Refs. [7, 12] in the framework of BBSvK counting. 
All these papers were written before the first unquenched lattice results appeared and hence 
do not use lattice data to fit their unknown free parameters. In both approaches the behavior 
of the scattering length was studied for a suitable range of the unknown parameters. Special 
attention was devoted to the extrapolations to the chiral limit. A more recent study can be 
found in Ref. [37] using the power counting of Ref. [18] and lattice data of the NPLQCD 
Collaboration. In the 1 So channel our results in the chiral limit indicate that the scattering 
length remains negative, thus the system is unbounded, coinciding with the predictions of 
mentioned previous works, albeit our value seems slightly smaller. In the 3 S\ channel our 
extrapolation of the scattering length to the chiral limit shows that it evolves from positive 
values at the physical pion mass to negative values, hence going from a bounded nucleon- 
nucleon system to an unbounded one. This is opposite to the results in Refs. [6, 36], and to 
those in Refs. [7, 12], for most of the parameter space, in which the scattering length remains 
positive in the whole range from the chiral limit to the physical pion mass. Nevertheless, in 
Refs. [7, 12] a behavior similar to the one we have obtained is observed in certain regions of 
the parameter space. In Ref. [37], the only one of the pevious works on the m 9 -dependence 
of the scattering lengths that has used lattice data, the 3 Si channel goes to negative values 
in the chiral limit, and overall presents a very similar result to ours. 
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6. Conclusions 



We have showed that certain classes of diagrams involving potential pion exchanges in loops 
with radiation pions can be summed up in the 3 Si channel. This is important because 
each of these exchanges introduces a parametric suppression of only 0(y/m n /rriN) that 
numerically turns out to be 0(1). The resummation is possible because after radiating 
a pion a nucleon-nucleon system in the 3 Si channel changes into the 1 So channel, and in 
this channel the one-pion exchange potential at high momentum transfer becomes a contact 
interaction. We showed that by performing dibaryon local field redefinitions we can get rid of 
the contact interaction, and hence the contribution of all diagrams involving these potential 
pion exchanges must be zero. We checked this cancellation by explicitly computing the 
diagrams and adding them up. 

Unfortunately, in the 1 Sq channel it has not been possible for us to compute the contri- 
bution of an arbitrary number of potential pions in a loop with a radiation pion. This is 
because after radiating a pion a nucleon-nucleon system in the 1 Sq channel changes into the 
3 Si channel, and in this channel the one-pion exchange potential at high momentum transfer 
does not reduce to a contact interaction anymore. However, similar arguments still apply 
to the diagrams with only one potential pion, which should then add up to zero. This is in 
contradiction with the results of Refs. [23, 25], and we have pointed out a possible source of 
the discrepancy in section 3.2. It is very likely that in the 1 S'o channel the perturbative series 
breaks down as in the 3 Si channel, which means that it is possible that our expressions for 
a So are missing large contributions, and hence, are unreliable. 

We have given chiral extrapolation formulas for l/a s ° and l/a Sl up to corrections of 
order O (m 3 /A^) depending on three independent free parameters. In section 5. we carried 
out a fit of these expressions to lattice data from the NPLQCD Collaboration [33, 34]. The 
results in Fig. 9, Table 3 and Table 4 show that our expressions for a iSl are much more 
compatible with lattice data than those for a s °, which could indicate that the missing, 
potentially large, contributions to a s ° previously mentioned do exist. Using these results to 
extrapolate the scattering lengths in the chiral limit, we obtain that a So keeps its negative 
sign, while a iSl changes from positive to negative. However, at this stage, lattice data 
sets available are rather small, with relatively large pion masses, and often computed using 
different approaches, making it difficult hold any strong statement in this respect. 
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Appendix 



The B((3i,(3 2 ) loop integrals used in section 3.2 are defined as follows: 



d d l q q 2 1 

(47r) d - 1 (q 2 + mj)ft (q 2 + ml - a 4 )^ ' 



(6.1) 
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